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FOREWORD 

This  is  the  Final  Technical  Report,  Part  I,  of  a  research  program 
sponsored  by  the  Aircraft  and  Crew  Systems  Technology  Directorate,  Naval 
Air  Development  Center,  Warminster,  Pennsylvania,  18974,  under  contract  No. 
N-62269-82-C-0250.  The  NADC  monitor  was  Mr.  Lee  W.  Gause.  Dr.  A.  S.  D.  Wang 
was  the  principal  investigator;  and  he  was  assisted  by  Miss  M.  Slomiana  and 
Mr.  Ronald  Bucinell,  Drexel  graduate  students. 

The  primary  objective  of  this  research  was  to  develop  a  computer  routine 
which  can  be  used  to  calculate  the  stress  field  of  a  laminate  containing  a 
propagating  delamination  crack.  By  incorporating  a  virtual  crack  opening 
procedure  in  the  routine,  the  crack-front  strain-energy  release  rate  can  be 
calculated  as  the  delamination  crack  is  incrementally  propagated.  Since  a 
general  delamination  crack  may  have  a  contoured  front,  this  routine  is  capable 
of  adjusting  to  new  boundary  configurations  and,  at  the  same  time,  yielding  an 
accurate  three-dimensional  stress  field  solution. 

The  timely  development  of  the  present  routine  represents  a  major  step 
toward  analytical  treatment  of  delamination  mechanisms,  perhaps  the  most 
damaging  failure  mode  in  structural  composite  laminates. 

In  Part  I,  an  overview  of  the  energy  release  rate  method  is  presented, 
as  it  is  applied  to  describe  delamination  growth  in  laminates  subjected  to  both 
static  and  fatigue  loads.  One-dimensional  growth  (free-edge  delamination)  as 
well  as  two-dimensional  growth  (coun toured  delamination)  problems  are  studied. 
All  results  are  correlated  with  past  and/or  present  experimental  data. 

Part  II  presents  details  of  the  computer  code  used  to  perform  these  cal¬ 
culations.  It  is  contained  in  the  report  NADC- 840 18- 60 . 
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INTRODUCTION 


The  delamination  problem  is  one  of  the  most  important  problems  affecting 
the  use  of  structural  composite  laminates.  It  is  often  encountered  as  the 
most  damaging  laminate  failure  mode.  The  mechanisms  of  delamination  are  ex¬ 
tremely  complex  rendering  a  satisfactory  analytical  treatment  practically 
un tractable. 

As  far  back  as  in  the  late  1960 's,  laboratory  coupon  tests  had  established 
that  the  tensile  strength  of  a  laminate  depends  on  the  laminate  stacking  se¬ 
quence  [1].  The  stacking  sequence  could  determine  whether  or  not  a  free  edge 
delamination  would  occur  under  the  applied  load.  If  a  free  edge  delamination 
occured,  it  would  rapidly  propagate  under  increasing  load,  causing  premature 
failure.  This  unique  phenomenon  can  be  present  not  only  in  laminates  under 
tensile  loads  [2],  but  also  in  laminates  under  compressive  loads  [3],  so  long 
as  the  laminate  stacking  sequence  is  such  that  free  edge  delamination  is  in¬ 
duced  by  the  applied  load. 

A  relatively  slower  delamination  growth  can  also  be  induced  by  cyclic 
loads,  be  it  tensile  and/or  compressive  [3].  Onset  of  delamination  can  occur 
at  cyclically  applied  stress  levels  much  below  the  threshold  stress  determined 
by  static  load.  Generally,  delamination  growth  considerably  shortens  the 
laminate's  fatigue  life. 

Regardless  of  the  nature  of  loading,  edge  delamination  has  been  thought  to 
be  driven  by  the  interlaminar  stresses  which  exist  near  the  laminate  edge 
boundaries.  Early  analytical  efforts  have  concentrated,  therefore,  on  the 
determination  of  the  free  edge  interlaminar  stresses.  These  are  highly  con¬ 
centrated  near  the  intersections  of  the  free  edge  and  the  lamina  interfaces 
[4-6].  In  fact,  these  intersection  points  are  singular  points  [7],  where  the 
stress  field  is  unbounded.  Because  of  mathematical  difficulties,  rigorous 
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free  edge  stress  analyses  have  been  limited  to  the  special  problem  of  a  long, 
uniform  width  and  symmetrically  stacked  laminate  subjected  to  uniaxial  tension. 

Nevertheless,  these  analyses  have  helped  to  unveil  the  complicated  nature 
oi  the  delamination  mechanisms.  In  particular,  knowledge  of  the  interlaminar 
stresses,  their  signs  and  magnitudes,  may  determine  whether  or  not  a  free 
edge  delamination  will  occur. 

Laboratory  experiments  have  shown  that  the  growth  behavior  of  free  edge 
delamination  depends  profoundly  on  the  absolute  thickness  of  the  material 
layers  of  the  laminate  [8,  9].  For  example,  the  onset  of  delamination  and 
the  subsequent  growth  behavior  under  tensile  loading  was  significantly  dif¬ 
ferent  between  the  [±45/0/90]g  and  [±452/02/902 ]g  laminates,  even  though  the 
calculated  interlaminar  stresses  are  identical  both  in  sign  and  in  magnitude 
[8]. 

Other  complexities  have  been  observed.  The  most  perplexing  are  the  phe¬ 
nomena  of  failure  modes  interactions.  For  Instance,  consider  the  [±45/0/90] 

S 

laminate  under  uniaxial  tension.  In  some  tests,  free  edge  delamination  oc- 
cured  first,  its  growth  inducing  multiple  transverse  cracks  in  the  90-layer. 
Conversely,  90-layer  transverse  cracks  occured  first  in  other  tests  which  then 
precipitate  localized  delaminations  [9].  The  growth  of  the  latter  is  generally 
multi-directional,  resulting  in  a  contoured  delamination  front. 

All  these  difficulties  have  complicated  the  development  of  a  satisfactory 
analytical  model.  Even  with  a  considerable  experimental  base  emcompassing  the 
various  facets  of  delamination,  no  rational  methodology  which  generally  describes 
the  onset  and  the  growth  of  delamination  had  been  devised. 

In  a  series  of  recent  papers  [10  -14],  Wang,  Crossman,  et.  al  successfully 
applied  an  energy  method  to  model  the  onset  and  growth  behavior  of  a  class  of 
simple  matrix  cracking  problems.  The  type  of  cracks  considered  included  the 

-2- 
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afore-mentioned  90-layer  transverse  cracks  and  free  edge  delamination  induced 
by  statically  applied  tensile  and  compressive  loads.  Their  model  proved  to  be 
generally  applicable  when  these  two  types  of  cracks  occurred  individually 
in  a  given  laminate. 

This  paper  extends  the  application  of  the  energy  method  to  fatigue  load 
induced  delamination  growth,  and  the  more  complex  contour  delamination  growth 
caused  by  two  interacting  flaws.  Correlations  between  the  analytical  and  the 
experimental  results  are  also  presented. 
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THE  ENERGY  RELEASE  RATE  METHOD 

The  energy  method  is  based  on  the  fundamental  concept  of  classical  frac¬ 
ture  mechanics.  Essentially,  it  involves  nothing  more  than  a  judicious  use 
of  the  crack  growth  criterion: 

G(tf,  a)  =  Gc  (1) 

Equation  (l)  is  the  well  known  Griffith  criterion  for  a  brittle  crack 
propagating  in  a  self-similar  manner  inside  a  homogeneous  solid  £15].  It 
actually  describes  the  energy  balance  at  incipient  propagation  for  a  crack 
having  a  length  "a"  and  under  the  applied  load  tf.  But  unlike  in  the  classical 
text-book  problems  where  the  crack  geometry( including  the  crack  length  a)  and 
the  material  fracture  toughness  Gc  are  well  defined  apriori ,  these  quantities 
are  either  generally  unknown  or  difficult  to  determine  for  the  types  of 
cracks  found  in  fibrous  composites.  In  addition,  mathematical  computation 
of  the  function  G(tf,a)  for  a  crack  propagating  in  multi-layered  anisotropic 
media  could  also  become  intractable. 

In  the  Wang-Crossman  approach,  these  limitations  were  largely  overcome 
through  the  use  of  the  following  fundamental  concepts  as  well  as  recent  de¬ 
velopments:  (a)  the  fundamental  postulate  of  the  "effective"  material  flaw 
distribution  which  enables  the  definition  of  the  initial  crack  length  "a" 
as  a  random  variable,  (b)  the  expanded  understanding  of  the  quantity  Gc  as 
a  local,  mode-dependent  and  directionally  variable  material  property,  (c)  the 
improvement  of  various  damage  detection  methods  which  provided  a  visual  des¬ 
cription  of  the  cracking  process  as  a  continous  function  of  loading  and/or 
time  under  load,  and  finally,  (d)  the  development  of  suitable  numerical  pro¬ 
cedures  to  compute  G(tf,a)  for  cracks  propagating  in  mid ti -layered  anisotropic 
solids.  A  brief  discussion  on  some  of  the  above  aspects  is  given  below. 
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The  Concept  of  Effective  Flaw  Distribution. 

In  laminate  stress  analysis,  the  fundamental  approach  has  been  the  so- 
called  ply-elasticity  theory  £l6].  The  theory  pertains  primarily  to  laminates 
that  are  made  of  unidirectional  plies  stacked  together,  with  each  ply  oriented 
in  some  designed  direction.  The  individual  plies,  though  consisting  fibers 
and  matrix  phases,  are  assumed  macro scopically  homogeneous  and  having  uni¬ 
form  material  properties.  This  assumption  enables  characterization  of  the  uni¬ 
directional  plies  by  two  sets  of  material  constants.  The  first  set  consists 
of  the  ply  deformation  constants,  while  the  second  set  is  associated  with  the 
ply  strength  constants.  For  example,  the  stiffness  matrix  in  the  gen¬ 
eralized  Hooke's  law  and  the  strength  constants  F.  .  and  F.  in  the  Quadratic 

*  J  J 

failure  criterion  are  among  the  most  widely  used  ply  properties  Cl6j: 


d.  =  G,  .e  . 
i  ij  J 

i» j  =  1.6 

(2a) 

F,  .d.d.  +  Fjd.  ~  1 
lj  i  J  J  J 

i.j  =  1,6 

(2b) 

It  is  noted  that  the  constants  F^,  etc.  are  regarded  as  uniform 

property  of  the  ply  material,  and  they  are  usually  determined  experimentally 
from  a  specimen  having  finite  volume.  In  particular,  failure  criterion  such  as 
expressed  in  Equation  (2b)  determines  the  first  and  only  failure  of  the  ply. 

As  has  been  discussed  earlier,  multiple  failures  in  a  certain  ply  can  occur 
at  different  stress  levels  when  that  ply  is  a  part  of  a  laminate.  Hence,  ply- 
strength  cannot  be  regarded  as  a  constant  material  property  in  this  context. 

The  conceptual  argument  for  multiple  ply  failures  in  laminates  is  that 
randomly  distributed  material  flaws  exist  in  the  micro structure  of  the  uni¬ 
directional  plies  and  in  the  ply-to-ply  interfaces.  The  true  physical  charac¬ 
teristics  of  the  flaws  are  believed  to  be  inherent  to  the  basic  ply  material 
itself,  and  possibly  also  to  the  laminate  fabrication  process,  post-process 
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handling,  etc.  Realistically,  the  true  identity  and  the  exact  physical  effect 
of  the  individual  flaws  cannot  be  determined  within  the  context  of  the  macro¬ 
theory  of  ply-elasticity.  One  useful  alternative,  however,  is  to  introduce 
the  concept  of  "effective"  material  flaws. 

This  concept  theorizes  that  a  characteristic  distribution  of  "effective" 
flaws  can  be  determined  for  a  given  representative  volume  of  the  basic  ply 
material;  each  of  the  effective  flaws  has  a  definite  size  as  well  as  location. 
Though  hypothetical  in  nature,  the  effective  flaws  must  be  so  characterized 
as  to  represent  the  aggregate  effect  of  local  microflaws  and  inperfections 
at  the  ply  level.  In  this  context,  each  effective  flaw  is  further  assumed  to 
act  like  a  small  crack;  and  it's  growth  behavior  may  be  descrlbable  by  the 
Griffith  criterion  (l).  Then,  each  of  the  effective  flaws  can  individually 
become  a  sublaminate  crack  at  some  critical  stress  level;  and  they  form  col¬ 
lectively  multiple  sublaminate  cracks  in  the  course  of  loading. 

Once  the  size,  the  location  and  the  geometrical  character  of  an  effective 
flaw  is  defined  in  relation  to  the  layering  structure  of  the  considered  la¬ 
minate,  the  flaw( crack)  tip  strain  energy  release  rate  G(d,a)  can  be  calcu¬ 
lated  by  some  suitable  computational  method. 

The  Crack-Closure  Procedure  for  G(d,a). 

For  a  line  crack  in  two-dimensional  elastic  solids,  the  stress  field  near 
the  crack  tip  is  singular  in  nature.  The  classical  methods  of  solution  have 
been  to  solve  for  the  near-field  stresses  by  means  of  a  boundary- value  problem 
in  the  theory  of  elasticity  £l?3.  These  stresses  are  generally  obtained  for 
three  distinct  modes  of  the  crack  actions  the  opening  mode-I ,  the  sliding  mode- 
II  and  the  anti -plane  mode-III.  For  each  mode  there  is  an  associated  stress 
intensity  factor,  Kj,  and  The  latter  are  functions  of  the  crack 

size  "a"  and  the  applied  far-fleld  stress  d  £18].  The  quantity  G(tf,a)  can  be 


-6- 
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related  directly  to  K^,  and  Kjjj*  This  procedure,  however,  can  become 
mathematically  intractable  for  problems  that  involve  complex  geometry,  such 
as  multiple  cracks  in  layered  composites. 

Other  solution  techniques  for  G(cf,a)  are  available,  notably  the  well-known 
J-integral  method  £l9,20]  and  the  method  of  crack-closure  by  Irwin  £21].  In 
particular,  Irwin's  method  recognizes  the  equality  between  the  elastic  strain 
energy  which  is  released  during  a  virtual  crack  extension,  and  the  work  done 
in  closing  it  again.  The  closing  of  the  virtual  crack  extension,  Aa,  yields 
the  solution  for  the  required  surface  traction  vector  5  over  Aa.  And,  the 
work  done  is 

,  Aa 

AW(d,a)  =  ~2~  f  (3.Au)  da  (3a) 

where  Au  is  the  relative  crack  surface  displacement  vector  over  Aa.  Note  that 
both  3  and  Au  depend  on  the  crack  size  "a’’  and  the  applied  far-field  stress  d. 
Then,  from  the  elastic  equivalence  the  crack  extension  rate  of  strain  energy 
release  can  be  expressed  by 

G(<*,a)  =  (AW/Aa)  (3b) 

The  advantage  of  expressing  G  in  the  form  of  (3a)  and  (3b)  is  that  the 
vector  product  (3-Au)  yields  a  sum  of  three  scalars  which  associate  directly 
with  the  three  modes  of  the  cracking  action.  Hence,  (3b)  gives  the  sum  of  G^ , 
Gjj  and  Gjjj*  In  addition,  the  integral  in  (3a)  can  easily  be  evaluated  nu¬ 
merically  by  discretization  of  the  solution  field  which  contains  the  considered 
crack.  For  example,  the  traction  vector  3  and  the  crack  extension  displacement 
vector  Au  may  be  approximated  by  the  nodal  forces  and  nodal  displacements, 
respectively,  in  a  finite  element  representation. 

Rybicki  and  Kanninen  £22]  first  suggested  the  finite  element  scheme  for 
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simulating  uniform  free  edge  delamination  in  a  laminate  subjected  to  uniaxial 
extension.  Wang  and  Crossman  QlO^  extended  the  scheme  in  a  generalized  plane 
strain  formulation  allowing  warping  of  the  laminate's  cross-section.  Their 
finite  element  routine  has  been  used  successfully  for  transverse  cracking  and 
free  edge  delamination  simulations  C 10-14], 


The  Physical  Charateri sties  of  G 

The  right-hand  side  of  Equation  (l)  is  a  material  property;  the  term  Gc 

represents  the  material  resistance  against  crack  propagation,  which  must  be 

measured  experimentally.  Generally  speaking,  the  measured  Gc  depends  on  the 

local  morphology  of  the  fractured  material  surface.  It  is  usually  regarded  as 

a  constant  in  homogeneous,  isotropic  materials.  For  fibrous  composites,  however, 

the  fractured  surface  morphology  can  differ  greatly  depending  on  the  local 

matrix/fiber  reinforcement  geometry  which  can  be  both  spatial  and  directional 

dependent.  But,  on  the  basis  of  ply-elasticity  from  which  the  fracture  growth 

criterion  (l)  is  applied,  such  local  surface  morphology  cannot  be  concurrently 

distinguished.  For  this  reason,  laboratory  measured  Gfi  in  matrlx-dominate 

crack  growth  has  been  found  to  depend  on  the  propagation  direction  and  the 

local  fiber  and/or  ply-interface  orientation  C23l* 

For  example,  Cullen  \_Zk~\  and  William  considered  two  different  cases 

of  mode-I  delamination  growth  in  unidirectional  laminates,  see  illustrations 

in  Figure  1.  The  first  case  concerned  the  crack  propagating  in  the  fiber 

direction,  while  the  second  case  concerned  the  crack  propagating  normal  to 

the  fiber  direction.  The  actual  fractured  surface  morphology  of  the  two  cases 

was  found  to  exhibit  a  considerable  amount  of  differences.  And,  as  a  result, 

the  respective  measured  G  also  differed  significantly.  Form  the  view  point 

c 

of  ply-elasticity,  however,  these  two  cases  should  not  display  any  difference. 

Another  aspect  of  the  measured  Gc  is  related  to  mixed-mode  cracks  that 
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contain  both  opening  and  shearing  modes.  Apparently,  the  total  strain  energy 
dissipated  in  a  mixed-mode  crack  extension  is  greater  than  in  a  purely  mode-I 
crack  extension.  The  increase  in  the  measured  Gc  is  generally  attributed  to 
crack- tip  matrix  yielding  under  shear  deformation.  In  several  recent  experi¬ 
ments  on  graphite-epoxy  composites,  the  total  G  for  mixed-mode  matrix  cracks 

c 

was  found  to  increase  with  Gjj/Gj  ratio  [[26-283 .  Consequently,  if  the  growth 
criterion  (l)  is  used  for  mixed-mode  cracking,  Gc  is  multiple-valued;  see, 
e.g.  the  curve  shown  in  Figure  2. 

Although  other  alternate  criteria  for  mixed-mode  cracks  have  been  used  in 
various  crack  growth  situations,  which  normally  require  more  than  one  constant 
(e.g.  Gjc,  etc.),  questions  surrounding  the  existence  of  and  the 

determination  method  for  the  constants  have  not  been  resolved  for  cracks 
propagating  in  composites.  Apparently,  the  exact  mechanisms  in  mixed-mode 
cracking  must  be  examined  carefully  at  the  microscopic  scale.  Such  a  topic 
still  remains  an  area  of  future  research. 

In  view  of  the  foregoing  discussions,  it  is  clear  that  the  classical 
criterion  as  expressed  in  Equation  (l)  cannot  be  applied  straight- forwardly 
for  sublaminate  crack  growth  problems  in  composite  laminates.  Judicious  inter¬ 
pretations  for  each  of  the  elements  in  the  equation  must  be  made  in  order  to 
develop  a  viable  analytical  model  for  the  cracking  phenomenon  as  observed  at 
the  macroscopic  scale. 

As  has  been  mentioned  previously,  the  energy  approach  has  been  applied 
successfully  to  describe  the  formation  process  of  multiple  transverse  cracks 
in  laminates  loaded  monotonically  or  cyclically  L293*  The  same  approach  also 
proved  applicable  to  model  uniform  free  edge  delamination  growth  when  it  is 
not  interacting  with  other  modes  of  cracking  C3()D» 

Some  key  features  of  the  free  edge  delamination  model  will  be  presented 
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next,  along  with  an  extended  model  for  fatigue  load  induced  delamination  growth. 
Both  one-dimensional  and  two-dimensional  growth  problems  will  be  analyzed  and 
compared  to  experimental  case  study  results. 
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ONE  -  DIMENSIONAL  GROWTH  PROBLEMS 

Consider  the  laminate  having  a  finite  width  and  subjected  to  uniform 
axial  tensile  strain  ex;  see  inset  in  Figure  J.  Let  the  laminate  stacking 
sequence  be  such  that  tensile  interlaminar  normal  stress  6^  is  induced  near 
the  laminate  free  edge  region.  Assume  that  an  effective  interlaminar  flaw 
having  a  random  size  "a"  exists  along  the  edge  as  shown.  When  the  free  edge 
stress  field  reaches  a  certain  critical  state,  the  effective  flaw  will  propa¬ 
gate  and  form  a  crack  opening  on  the  free  edge.  This  event  then  constitutes 
the  "onset"  of  free  edge  delamination,  and  the  concurrent  stress  state  is 
termed  the  threshold  stress. 

For  clarity,  let  it  be  assumed  that  the  particular  interface  in  which  de¬ 
lamination  will  occur  is  known  apriori .  Then,  the  flaw( crack)  tip  energy  re¬ 
lease  rate  G(ex,a)  for  the  assumed  effective  flaw  "a”  can  be  calculated  by 
the  finite  element  crack-closure  procedure.  For  a  given  value  of  e^,  G  will 
vary  with  "a".  A  typical  G-curve  is  shown  in  Figure  3.  It  is  seen  that  G  rises 
sharply  from  zero  at  a  =  0,  and  reaches  an  asymptotic  value  G^  at  a  £  a^ 

It  is  of  interest  to  note  that  although  a  stress  singularity  may  exist 
at  the  free  edge/ply  interface  intersection  £7],  a  non- trivial  G  could  be 
determined  only  when  an  effective  flaw  of  size  a  >  0  is  introduced  near  the 
ply  interface. 

As  for  the  value  of  am,  it  is  usually  of  the  order  of  the  ply  thickness. 
Hence,  the  value  of  G^  can  be  affected  by  the  ply-thickness  parameter  £l4]. 

With  linear  elastic  plies,  it  is  often  convenient  to  express  G  in  ex¬ 
plicit  terms  of  the  applied  strain  ex»  For  instance 

G(ex,a)  =  Ce(a) -t-e^  (4) 

where  C  is  a  function  of  the  flaw  size  (a/t)  only;  and  t  is  a  chosen  length 
© 

parameter,  usually  taken  as  one  ply  thickness. 
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Similarly,  if  the  laminate  is  subjected  to  a  uniform  temperature  load,  AT, 
and  the  thermally  induced  free  edge  stresses  can  also  effect  delamination,  then 
the  associated  flaw  tip  energy  release  rate  G  is  expressible  as, 

G(AT,a)  =  CT(a).fdT2  (5) 

Where  C^,  is  also  a  function  of  the  flaw  size  (a/t)  only,  and  whose  functional 
shape  is  similar  to  that  shown  for  C0. 

Normally,  AT  is  determined  experimentally  to  represent  the  temperature  drop 
from  the  stress-free  tempera ture( :£  laminate  curing  temperature)  to  the  ambient 
service  temperature. 

The  mechanical  and  the  temperature  load  effects  on  the  laminate  stress 
field  are  linearly  additive j  and  the  total  energy  release  rate  for  the  combined 

loading  case  can  be  expressed  as, 

G(Sx,AT,a)  =  Uc~e  -ix  ♦  -AT]2-!  (6) 

If  the  size  of  the  edge  flaw  "a"  is  given,  then  the  onset  of  it's  propaga¬ 
tion  also  delamination)  is  defined  when  the  applied  laminate  strain  is  such 
that 

G(Sx,AT,a)  =  Gc  (7) 

With  "a"  being  random  valued,  however,  Equation  (7)  can  be  plotted  in  the 
ex  versus  "a"  plane.  Then,  given  a  range  of  possible  values  for  the  edge  flaw 
size,  a  corresponding  range  for  the  onset  load(i.e.  ex)  can  be  determined 
from  the  plot  as  shown  in  Figure  4.  In  particular,  if  for  all  probabilities 
the  size  of  the  edge  flaw  is  of  the  order  of  am  or  larger,  then  the  predicted 
load  range  for  onset  of  delamination  is  very  narrow;  in  fact,  it  becomes  prac¬ 
tically  a  constant,  see  Figure  4.  The  latter  represents  the  minimum  possible 
load  for  onset  of  delamination. 
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Illustrative  Example  -  Laminates  Under  Static  Tension. 

In  order  to  illustrate  the  use  of  the  eneigy  model  outline  above,  consider 

the  laminate  family  £±25/90,,]  ,  n  =  1,2  and  3»  The  laminates  were  made  using 

n  s 

AS- 3501-06  graphite-epoxy  system,  and  were  tested  at  room  temperature  under 
static  uniaxial  tension  (for  details,  see  Ref.  £30]).  Free  edge  delamination 
growth  was  induced  in  all  laminates  during  the  course  of  loading.  Figure  5 
shows  the  experimental  plot  for  delamination  size  versus  applied  laminate 
stress  In  each  of  the  three  types  of  laminates,  three  replicates  were  used; 
and  an  averaged  onset  load  for  delamination  could  be  determined  by  extrapola¬ 
tion  of  the  delamination  growth  curve. 

Post- test  examination  of  the  specimens  showed  that  delamination  in  the 

£+25/90]  laminate  was  contained  in  the  mid-plane,  while  delamination  in  the 
s 

£±25/90 „]  and  the  £±25/90 laminates  was  in  the  -25/90  interface.  Actually, 

C,  S  J)  S 

90-layer  transverse  cracks  in  the  latter  two  laminates  were  observed  to  precede 

free  edge  delamination.  When  delamination  occurred,  it  was  mostly  localized 

near  clusters  of  transverse  cracks  £30].  The  latter  is,  of  course,  a  case  of 

"crack  modes  interaction",  which  will  be  treated  later  in  this  paper  as  a  two- 

dimensional  contour  delamination  growth  problem. 

Only  in  the  £±25/90 ]  laminate  was  onset  of  edge  delamination  not  compli- 

s 

cated  by  transverse  cracks.  The  delamination  grew  uniformly  along  the  length 
of  the  specimen.  Hence,  the  fracture  model  for  one-dimensional  growth  as  out¬ 
lined  above  could  be  applied.  Using  the  base  properties  of  the  unidirectional 
ply  (see  Ref.  £30]  for  values)  in  the  energy  release  rate  calculation,  the 
coefficient  functions  appearing  in  Equation  (6)  were  generated  for  an  assumed 
delamination  in  the  laminate  mid-plane  and  in  the  -25/90  interface.  These 
calculated  functional  curves  are  displayed  in  Figures  6  and  7,  respectively. 

It  is  noted  that  mid-plane  delamination  is  of  mode-I  due  to  symmetry  in 
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the  lamination  stacking  sequence,  while  -25/90  interface  delamination  is  of 
mixed-mode  (I  and  III).  In  the  present  example,  G^^  is  about  0.8  times  of  Gj. 

From  Figures  6  and  7,  it  is  seen  that  the  total  energy  release  rate  for 
mid-plane  delamination  is  larger  than  that  for  -25/90  interface  delamination. 
Furthermore,  the  material  resistance  Gc  for  mode-I  growth  is  smaller  than  that 
for  mixed-mode  growth  (see,  e.g.  Figure  2).  Thus,  the  model  would  predict  a 
mid-plane,  mode-I  delamination  growth.  Using  the  maximum  values  of  Ce  and  CT 

in  Figures  6  and  7i  the  minimum  onset  load  (in  terms  of  far- field  laminate 
strain  ex)  can  be  calculated  from  Equations  (6)  and  (7)l  thus, 

<5*>  min  =  0-53* 

In  computing  above,  Gj  =  175  J/m2  (1.0  in-lb/in2),  t  =  0.132  mm  (0.005") 
and  AT  =  125°C  (225°F)  were  used,  see  Ref.  £30]. 

For  the  laminate  considered,  the  axial  modulus  Ix  =  64  Gpa  (9.3  msi) ;  then 
the  calculated  minimum  onset  far- field  stress  is, 

(3X)  mm  =  338  Mpa  (49.0  ksi) . 

Extrapolating  the  experimental  value  for  (3  )  from  Figure  5  yields  a 

x  cr 

range  from  331  to  345  Mpa  (48  to  50  ksi).  The  prediction  falls  in  this  range. 

Illustrative  Example  -  Laminates  Under  Static  Compression. 

To  further  assess  the  usefulness  of  the  energy  model,  consider  the  follow¬ 
ing  laminates  which  are  geometrically  similar »  £902/02/t452]s,  £o2/902/t452Dg 

and  C(0/90)2/(t45)2]g.  These  laminates  were  also  made  using  the  AS-3501-06 

graphite-epoxy  system;  but  they  were  tested  in  room  temperature  under  static 

compressive  loading.  In  each  test  case,  three  replicate  specimens  were  used. 

The  stacking  sequence  of  all  the  three  types  of  laminates  was  designed 

to  yield  a  tensile  interlaminar  normal  stress  <?z  near  the  free  edge  under  the 
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applied  compressive  load.  And,  indeed,  edge  delamination  in  the  laminate 
mid-plane  was  induced  in  all  test  cases  (detailed  test  data  may  be  found  in 
Ref.  t3  ]  Figure  8  shows  the  delamination  growth  as  plotted  against  the 

applied  stress  5^  for  each  of  the  three  types  of  laminates.  The  growth  curves 
indicate  that  onset  of  edge  delamination  was  followed  by  a  rapid  growth,  which 
caused  the  final  failure  of  the  laminates.  Note,  in  particular,  that  the  onset 
load  for  delamination  differed  appreciably  among  the  three  types  of  laminates, 
though  they  are  in  many  ways  similar  to  one  another. 

A  free  edge  stress  analysis  was  performed  in  order  to  discern  any  major 

differences  which  may  be  present  in  their  free  edge  stress  distributions. 

Figures  9»  10  and  11  show  the  through-thickness  d  -distribution  for,  respec- 

z 

tively,  C902/02/±452]s,  C02/902/±452]s  and  C(0/90)2/(±45)2]g.  Similar  re¬ 
distributions  are  shown,  respectively,  in  Figures  12,  13  and  14.  Since  d  and 

z 

T^z  are  predominant  on  both  the  mid-plane  and  the  45/-45  interface  (in  fact, 
both  stresses  are  singular  on  the  45/-45  interface) ,  it  is  not  immediately 
clear  in  which  plane  edge  delamination  would  occur. 

Hence,  an  energy  release  rate  calculation  was  subsequently  followed.  The 
energy  release  rate  coefficient  Cg  for  delamination  in  the  mid-plane  and  in 
the  45/-45  interface  was  computed  separately;  and  they  are  shown,  respectively, 
in  Figures  15  and  16.  The  corresponding  C^,  curves  were  not  shown  because  of 
their  relatively  small  values. 

Flom  the  (^-curves  shown  in  Figures  15  and  16,  it  may  now  be  inferred  that 
a  mode-I,  mid-plane  delamination  would  occur  in  all  the  three  types  of  laminates. 
This  is  so  because  the  computed  Ge~curve  for  mid-plane  delamination  is  higher 
than  the  C  -curve  for  delamination  in  the  45/-45  interface  plane. 

To  predict  the  actual  onset  load,  the  Cg-curve  for  mid-plane  delamination 
must  be  used.  Thus,  as  was  done  in  the  previous  example,  the  minimum  onset 
strain  e^  can  be  calculated  using  Equations  (6)  and  (7).  The  results  are, 
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Laminate  Type 
C902/02/*452]s 
Co2/9o2/i45238 
C(0/90)2/(±45)2]s 


Minimum  Onset  <3 
— - x 

-303  Mpa  (-44  ksi) 
-352  Mpa  (-31  ksi) 
-40?  Mpa  (-59  ksi) 


The  corresponding  experimental  values  may  be  determined  from  Figure  8, 
yielding  the  following  respective  load  ranges:  290  -  310  Mpa*,  331  -  358  Mpa? 
and  400  -  427  Mpa.  The  predicted  results  are,  again,  conforming  with  the  ex¬ 
perimental  results. 

In  the  example  problems  discussed  in  this  section,  the  question  regarding 

the  size  "a"  of  the  assumed  effective  flaw  was  not  dealt  with  directly.  The 

assumption  was  simply  that  along  the  free  edge  of  the  laminate,  the  worst  flaw 

has  a  size  "a"  which  is  at  least  equal  or  greater  than  a  ;  and  a  is  of  the 

mm 

order  of  the  layer  thickness  in  which  the  assumed  delamination  is  contained. 
The  predJ  ’ted  onset  load  for  delaminatlon  is  thus  the  minimum  possible  load. 
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ONE  -  DIMENSIONAL  GROWTH  UNDER  FATIGUE  LOAD 

It  is  well  known  that  cyclic  fatigue  loads  also  induce  sublaminate  matrix 
cracks,  including  many  forms  of  delaminatlon  growth.  In  fact,  such  cracks  can 
be  induced  and  propagated  at  fatigue  stress  levels  much  lower  than  the  threshold 
stress  determined  statically.  Generally,  fatigue  induced  sublaminate  cracks 
axe  similar  to  that  induced  by  static  load,  although  the  respective  mechanisms 
may  be  fundamentally  different. 

In  a  recent  paper  by  Wang,  et.  al.  ^293,  the  formation  process  of  transverse 

cracks  in  £0/90 ~\  type  laminates  subjected  to  constant  amplitude  fatigue  loading 
s 

was  investigated.  It  was  demonstrated  that  given  a  fatigue  load  amplitude,  say 
a  characteristic  crack-density  versus  fatigue-time  (cycles)  relationship 
could  be  obtained  experimentally,  as  well  as  analytically  by  means  of  a  sto«- 
chastic  simulation  model.  Here,  Wang,  et.  al.  invoked  once  again  the  concept 
of  effective  flaw  distribution.  It  was  argued  that  inherent  material  flaws  do 
propagate  and  coalesce  under  cyclic  stress j  however,  the  aggregate  effect  within 
a  representative  macro-volume  may  be  described  by  the  propagation  of  an  assumed 
effective  flaw  having  the  size  "a".  In  this  context,  Wang,  et.  al.  assumed 
that  the  driving  force  that  propagates  "a"  during  each  cycle  of  loading  is  re¬ 
lated  to  the  instantaneous  flaw  tip  energy  release  rate  G(<J^,a),  which  can  be 
calculated  as  in  the  case  of  statically  loaded  problems.  Specifically,  a  simple 
power  law  similar  to  that  used  for  fine-grained  metals  C3lD  was  used  to  describe 
the  growth  rate: 


=  a  [G(tff,a)/Gc]  P  (8) 

where  N  is  the  fatigue  cycle  and  a  and  p  are  two  disposable  constants  to  be 
determined  empirically. 

It  is  noted  that  the  quantity  G(tf^,a)  in  Equation  (8)  already  includes 
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the  laminate  stacking  sequence,  lamina  thicknesses,  laminate  shape  features, 
the  individual  flaw  geometry,  the  nature  of  the  applied  load,  etc.  And,  fur¬ 
thermore,  the  quantity  G/Gc  represents  the  crack-driving  force  relative  to  the 
material's  resistance.  Hence,  the  fatigue  constants  a  and  p  are  expected  to  be 
material  dependent  only.  In  particular,  the  exponent  p  is  the  most  important 
parameter  for  the  crack  growth  rate. 

At  this  point,  it  seems  appropriate  to  review  briefly  some  related  works 
concerning  the  growth  rate  equation  (8)  and  the  constants  associated  with  that 
equation. 

Wilkins  C23l  conducted  fatigue  delamination  growth  experiments  using  split¬ 
ting  cantlliver  beams,  and  measured  the  growth  rate  as  a  function  of  the  crack 
tip  energy  release  rate  G(df,a).  Their  test  results  support  the  power  law  for 
delamination  growth,  such  as  expressed  in  Equation  (8).  By  taking  the  logarithm 
on  both  sides  of  (8),  a  linear  relation  between  da/dN  and  G(tff,a)  is  obtained; 
and  the  value  of  p  is  the  slope  of  the  straight  line.  Figure  1 7  shows  one  such 
experimentally  measured  da/dN  versus  G(tff,a)  plot  from  a  mode-I  delamination 
growth  test.  Figure  18  shows  one  for  a  mixed-mode  delamination  growth  case. 
Wilkins'  results  generally  confirm  the  adequacy  of  the  assumed  growth  rate  law 
for  both  mode-I  and  mixed-mode  cases.  However,  the  exponent  p  for  mode-I  growth 
is  found  to  be  two  to  three  times  larger  than  that  for  the  mixed-mode  case. 
Since  the  value  of  G/G  in  (8)  is  normally  smaller  than  unity,  a  larger  and 

C 

positive  p  implies  slower  growth  rate.  Wilkins  £23]  suggested  that  mode-I  de- 
lamination  may  be  predominantly  a  static  load  failure,  while  mixed-mode  de- 
lamination  may  be  predominantly  a  fatigue  load  failure. 

In  an  another  study  by  Odom  and  Adams  C32],  the  fatigue  property  of  the 
3501-06  epoxy  resins  was  investigated.  Two  types  of  fatigue  tests  were  per¬ 
formed!  a  tension- tension  (R  =  0.1)  fatigue  using  a  flat  dog-bone  coupon,  and 
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a  torsion- torsion  fatigue  using  a  round  dog-bone  coupon  (also  for  R  =  0.1). 

Both  tests  were  conducted  under  ambient  room  temperature  condition,  and  all 
tests  were  carried  to  the  final  failure  of  the  specimen.  The  conventional  S  -  N 
curves  were  obtained. 

Let  it  be  assumed  that  some  dominant  flaw  of  size  "a"  existed  on  the  sur¬ 
face  of  the  specimen,  and  it  propagated  stably  under  the  cyclic  stressing.  The 
growth  rate  of  the  propagating  flaw  under  the  applied  fatigue  amplitude  is 
governed  by  Equation  (8).  Since  the  energy  release  rate  at  the  flaw  tip  can 
be  expressed  in  terms  of  Equation  (8)  may  be  written  in  the  form 

ff  =  SUM^)2*  (9) 

For  constant  cf^.,  a  simple  Integration  of  (9)  gives, 

F(a)  =  N-(df)2p  (10) 

where  F(a)  is  a  known  function  of  the  flaw  size  "a".  Assume  that  final  failure 
occurred  when  a  acr  at  N  =  Ncr*  The  latter  was  measured  experimentally  and 
F(a  )  became  some  constant.  If  the  growth  rate  law  (9)  is  true,  then  a  linear 
relation  between  log(d-)  and  log(N  )  is  obtained;  the  slope  of  the  straight 

X  cx 

line  yields  the  value  of  2p. 

Figure  19  shows  the  log(tf-)  versus  log(N  )  plots  for  the  tension  and  the 

x  cr 

torsion  fatigue  tests.  These  results,  again,  support  the  assumed  growth  rate 
law. 

Note  that  failure  in  the  tension  fatigue  case  is  of  mode-I  in  nature.  The 
associated  fatigue  exponent  p  is  about  14,  see  Figure  19.  On  the  other  hand, 
failure  in  the  torsional  case  is  due  mostly  to  shearing  and  some  amount  of 
axial  tension  due  to  a  small  pre-tension  applied  to  the  test  specimen  (see, 

Ref.  C32D) •  The  latter  may  be  considered  as  a  mixed-mode  failure  case.  The 
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associated  value  for  p  is  about  7*5«  Although  these  values  are  found  for  the 
pure  epoxy  resin,  they  are  consistent  with  those  found  for  composites  made  of 
similar  resins.  Wilkins  C23I .  for  example,  found  the  value  of  p  for  mode-I 
growth  a  range  from  20  to  30i  and  for  mixed-mode  growth  a  range  of  7  to  8. 

Wang,  et.  al.  C29D  used  p  =  20  in  their  mode-I  transverse  cracking  simulation. 

Both  used  the  same  AS-fiber  and  the  3501-06  resin  system. 

Thus,  it  appears  that  the  delami nation  growth  rate  equation  (8)  can  be 
safely  used  to  simulate  it's  growth  behavior,  although  the  exact  physical 
mechanisms  involved  in  the  process  could  not  be  articulated.  The  latter  remains 
outside  the  scope  of  the  present  paper. 


Illustrative  Example  -  Laminates  Under  Tension  Fatigue. 

As  an  illustrative  example,  consider  the  C±25/90]s  laminate  discussed  earlier 
in  the  static  loading  case.  Recall  that  the  static  threshold  stress  for  delamina¬ 
tion  was  about  338  Mpa  (49  ksi)  and  the  growth  was  basically  mode-I  in  nature. 

As  the  laminate  was  fatigue  loaded  at  stress  amplitudes  below  the  static  threshold 
stress,  free  edge  delamination  in  the  -25/90  interface  was  observed.  This,  of 
course,  is  a  mixed-mode  growth  case. 

Although  the  computed  crack  tip  energy  release  rate  G(df,a)  is  larger  for 
mode-I  growth  than  for  mixed -mode  growth,  the  latter  is  nevertheless  predominant 
under  cyclic  stressing.  This  finding  conforms  with  Wilkins'  contention  that 
mixed-mode  delamination  growth  is  usually  a  fatigue  related  failure  mode. 

For  the  considered  laminate,  the  energy  release  rate  coefficients  Cg  and  CT 
associated  with  delamination  in  the  -25/90  interface  are  shown  in  Figures  6  and 
7.  Then,  by  means  of  the  growth  rate  equation  (8),  the  delamination  growth  can 
be  predicted  in  terms  of  the  fatigue  amplitude  d^  and  the  cycle  number  N. 

Figure  20  shows  the  experimental  delamination  growth  curves  under  fatigue 
amplitudes  d^.  *  207  Mpa  (29  ksi)  and  276  Mpa  (40  ksi).  The  solid  lines  sure  the 
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2  2 

corresponding  predictions  calculated  based  on  Gc  =  228  j/m  (1.3  in-lb/in  ) 
for  mixed-mode  growth,  and  the  values  for  a  and  p  were  taken  as 

a  =  1.016  mm  (0.04  inch);  p  =  8.3  (ll) 

It  is  seen  that  the  predictions  are  satisfactory.  Note  also  that  the 
value  for  p  used  here  falls  within  the  range  found  by  Wilkins  Q23D  and  by 
Odom,  et.  al.  £32^. 

Illustrative  Examples  -  Laminates  Under  Compression  Fatigue. 

Consider  further  as  examples  the  compressively  loaded  laminate  family 
[902/02A452]s.  C02/902/±452]s  and  £(o/90)2/(±45)23s.  The  respective  static 
threshold  stresses  for  edge  delamination  were  determined  previously  at  -303 
Mpa  (-44  ksi) ,  -352  Mpa  (-51  ksi)  and  -40?  Mpa  (-59  ksi) ;  see  Figure  8.  When 
these  laminates  were  tested  under  compression  -  compression  fatigue  loading 
(R  =  0.1),  stable  growth  of  edge  delamination  was  recorded  at  various  cycle 
intervals  (for  details,  see  Ref.  £3]) •  The  delamination  crack,  however,  occurred 
mostly  in  the  -45/45  interface,  indicating  a  mixed-mode  growth  case.  Tills  find¬ 
ing  contrasts  once  again  with  the  results  of  the  static  loading  case  where  the 
edge  delamination  was  essentially  of  mode-I  (and  confined  in  the  laminate's 
mid-plane) . 

To  simulate  the  mixed -mode  growth,  the  energy  release  rate  curves  for  edge 
delamination  in  the  -45/45  interface  were  used  in  conjunction  with  the  growth 
rate  equation  (8).  Figure  21  shows  a  comparison  between  the  experimental  and 
the  predicted  results  for  the  £902A>2/±452]S  and  £c>2/902/+452]s  laminates, 
both  under  <3^  =  -221  Mpa  (-32  ksi).  The  predicted  results,  shown  in  solid  lines, 
were  calculated  using  the  same  a  and  p  values  as  listed  in  Equation  (ll).  It 
is  noted  that  the  £902/o2/±452ls  laminate  was  the  weaker  under  static  loading, 
and  it  remained  the  weaker  under  the  fatigue  loading. 
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Similarly,  Figure  22  shows  a  comparison  between  the  Co?/90  /±45_]  and 
the  L( 0/90) 2/( ±45) 2]s  laminates  when  both  were  loaded  under  =  -2?8  Mpa 
(-40  ksi) .  Note  that  the  latter  laminate  displayed  a  distinctly  stronger  re¬ 
sistance  against  delamination,  as  it  also  did  when  under  statically  applied 
load. 

All  these  physical  features  pertaining  to  the  onset  and  growth  of  delamina 
tion  could  be  predicted  by  the  energy  model.  This  is  particularly  remarkable 
considering  the  scattering  nature  of  the  fatigue  results.  More  comparative 
studies  using  other  laminate  geometries  and/or  other  loading  combinations  have 
been  reported  in  Ref.  £30]  and  [33]. 
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TWO  -  DIMENSIONAL  GROWTH  PROBLEMS 

In  the  preceding  sections,  the  energy  release  rate  method  has  been  applied 
to  simulate  a  number  of  free  edge  delamination  growth  problems.  The  laminates 
studied  were  so  designed  and  loaded  as  to  induce  delamination  growth  without 
significant  interaction  with  other  cracking  modes.  Furthermore,  the  delamina¬ 
tion  was  confined  to  grow  uniformly  along  the  straight  edge  of  the  test  coupon, 
so  it  could  be  assumed  as  one-dimensional,  self-similar  growth.  Under  this 
assumption,  the  kinematics  of  the  crack  propagation  is  much  simplified,  with 
the  crack  front  being  represented  by  a  point  known  as  the  crack- tip,  and  the 
crack  size  having  the  magnitude  "a".  This  simplification  has  made  a  two-dimen¬ 
sional  stress  field  solution  possible,  from  which  the  crack-tip  energy  release 
rate  function  G(5,a)  could  be  calculated  by  a  numerical  technique. 

Generally  speaking,  one-dimensional  delamination  growth  is  a  very  special 
case,  and  it  happens  only  under  ideally  controlled  conditions.  Most  of  the 
problems  encountered  in  practice  are  localized  events  usually  caused  by  local 
defects.  The  growth  of  localized  delamination  is  almost  always  multi -directional 
in  nature. 

For  example,  consider  a  laminate  having  a  small  through-hole.  Upon  loading, 
localized  delamination  may  be  induced  near  the  curvilinear  edge  of  the  hole. 

In  this  case,  the  crack  front  will  be  some  line-contour,  and  the  associated 
growth  will  be  two-dimensional.  The  instantaneous  crack  size  "a"  is  now  the 
delaminated  area,  which  will  have  both  a  magnitude  and  a  definite  shape  contour. 

Clearly,  a  three  dimensional  stress  field  solution  is  needed  in  order  to 
calculate  the  crack  front  energy  release  rate  G(5,a) ,  with  "a"  being  some 
definite  line  contour. 

In  a  separate  effort,  Wang,  Kishore  and  Li  [.351  developed  a  3~D  finite 
element  routine  in  which  Irwin's  crack-closure  concept  was  incorporated  to 
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calculate  the  crack  front  energy  release  rate  along  a  prescribed  delamination 
contour.  Utilizing  this  3-D  finite  element  routine,  several  two-dimensional 
crack  growth  problems  have  been  simulated.  The  following  section  examines  one 
such  problem. 

Illustrative  Example  -  Growth  From  A  Pre-Planted  Delamination. 

Earlier  in  this  paper,  compression  induced  free  edge  delamination  in  three 
similar  quasi -iso tropic  laminates  was  examined  and  modeled  as  one-dimensional 
growth  problems,  see  Figure  8.  One  of  the  laminates,  Co2/9°2^ '-^2} s’  Was  later 
tested  with  an  implanted  delamination  at  it’s  mid-plane ,  see  insert  in  Figure 
23.  The  experimental  details  regarding  the  implantation  and  test  procedures 
were  reported  in  Ref.  C3U *  Figure  23  shows  the  test  results  from  three  specimens 
with  the  implant.  When  compared  with  specimens  without  the  implantation,  the 
laminates  with  the  implant  show  a  strikingly  different  growth  behavior.  Speci¬ 
fically,  the  laminates  with  the  implant  underwent  a  contoured  (2-D)  delamination 
growth,  while  the  ones  without  an  implant  experienced  more  or  less  uniform  (1-D) 
free  edge  delamination.  Although  the  contoured  growth  was  much  more  slow  and 
stable,  it's  early  initiation  had  actually  led  to  early  laminate  failure,  see 
Figure  23. 

In  order  to  predict  the  onset  of  the  contoured  delamination,  the  3-D  finite 
element  simulation  routine  developed  in  Ref.  E35D  was  used. 

But  before  preceding  with  the  growth  simulation,  a  stress  analysis  was 
performed  with  emphasis  placed  at  the  location  where  the  most  severe  stress 
concentration  is  present.  This  location  is  the  interface  corner  where  the  im¬ 
planted  delamination  front  Intersects  the  free  edge  of  the  test  coupon,  see 
Figure  24,  labeled  "interface  corner". 

It  is  recalled  that  the  laminates  without  the  implant  developed  a  tensile 
interlaminar  d^  stress  near  the  laminate  free  edge.  And,  for  a  far- field  strain 
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ex  -  -10”^  the  value  of  d^  near  the  free  edge  (acting  on  the  mid-plane)  is 
about  14  kpa  (2  psi) .  On  the  other  hand,  the  laminates  with  the  implant  deve¬ 
loped  a  strikingly  different  ^-distribution  in  the  same  laminate  mid-plane, 
see  Figure  24.  The  stress  near  the  "interface  corner”  is  in  fact  singular  in 

nature j  the  actual  value  of  d  near  the  vicinity  is  several  times  larger. 

z 

In  fact,  d is  not  the  only  stress  which  is  amplified  at  this  point  (the 

two  interlaminar  shear  stresses  are  zero  due  to  lamination  symmetry).  Figure  25 

shows  the  in-plane  stresses  d  ,  d  ,  and  T  along  with  d  that  are  acting  in 
*  x  y  xy  z  D 

the  45-layer  adjacent  to  the  laminate  mid-plane.  It  is  seen  that  all  of  the 
stresses  are  severely  concentrated  near  the  "interface  corner"  point.  These 
stresses,  particularly  d^,  are  obviously  responsible  for  the  observed  early 
growth  of  contoured  delamination. 

Returning  to  predicting  the  associated  delamination  growth,  one  must  begin 
with  the  calculation  of  the  crack  front  strain  energy  release  rate.  Since  the 
expected  delamination  will  start  from  the  "interface  corner”  point,  and  grow 
with  a  certain  contoured  front,  a  reasonable  sequence  of  nodal  release  may  be 
prescribed.  For  example,  Figure  26  (a)  shows  one  possible  nodal  release  sequence 
which  simulates  the  delamination  growth  process  (the  node  numbers  represent 
the  order  of  nodal  release) .  The  computed  energy  release  rate  coefficient  Cg 
as  a  function  of  the  correspondingly  delaminated  area  is  shown  in  Figure  26  (b). 

It  is  seen  that  the  delamination  growth  is  most  energetic  along  the  implant 
edge  as  compared  to  the  energy  released  along  the  laminate  free  edge.  As  the 
delaminated  area  becomes  larger,  the  available  energy  release  rate  becomes 
smaller.  This  implies  that  the  growth  process  is  essentially  stable,  as  it 
has  been  observed  in  the  experiment,  see  Figure  23. 

To  determine  the  onset  load,  the  maximum  value  of  Cg  in  Figure  26  (b)  will 
be  used  in  order  to  predict  the  minimum  onset  load.  By  disregarding  the  thermal 
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effect  on  delamination  (because  C^,  is  relatively  small),  the  maximum  strain 
energy  release  rate  can  be  expressed  by, 

2 

G  =  (C  )  (5  )  t  (12) 

max  e  max  x 

where  the  value  of  (Ge)max  corresponds  to  the  release  at  node  6,  see  Figure  26. 

Using  G  =  1?5  J/m2  (1,0  in-lb/in2)  for  mode-I  delamination,  and  t  =  0.13 
c 

mm  (0.005  in.),  the  calculated  minimum  onset  load  (the  applied  far- field  strain 
ex)  is  approximately, 

-6 

(®x>mln  *  28*>  1110  03) 

For  the  laminate  considered,  the  axial  stiffnes  is  E^  =  50.1  Gpa  (7.27  msi) 5 
it  follows  that  the  predicted  (minimum)  onset  fax-field  laminate  stress  is 

(3x)Bln  =  145  Mpa  (21’°  ksi)  (1*0 

This  predicted  value  is  just  about  the  same  as  the  experimental  value  which 
can  be  extracted  from  Figure  23. 
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CONCLUDING  REMARKS 

In  this  paper,  the  energy  release  rate  concept  was  applied  to  describe 
the  initiation  and  growth  processes  in  some  of  the  most  prevalent  matrix  cracks 
found  in  composite  laminates.  Several  experimental  case  studies  were  used  in 
order  to  illustrate  the  analytical  models.  At  the  conclusion  of  this  paper, 
the  following  remarks  are  in  order: 

1.  For  matrix  cracks  such  as  delamination,  the  energy  release  rate  method 
can  be  applied  within  the  frame  work  of  ply  elasticity  and  the  classical 
fracture  mechanics.  To  do  so,  all  variables  associated  with  the  analytical 
model  must  be  defined  or  determined  at  the  material  ply  level. 

2.  For  self-generated  matrix  cracks  (cracks  caused  by  inherent  material 
flaws) ,  a  concept  of  "effective**  flaw  can  be  used  in  place  of  the  conven¬ 
tional  strength  concept.  The  new  concept  is  consistent  with  the  fundamental 
postulations  in  the  classical  fracture  mechanics.  And,  by  defining  it  as 
an  "effective"  basic  material  property,  the  concept  is  also  consistent  with 
the  assumptions  of  ply  elasticity  (that  all  material  properties  are  defined 
at  the  ply  level) .  This  concept  provides  the  necessary  arguement  for  multi¬ 
ple  crack  formations  observed  in  laminates,  for  which  the  conventional 
constant  ply  strength  concept  cannot. 

3.  One-dimensional  delamination  growth,  such  as  uniformly  delaminated  free 
edge,  is  a  special  case.  Most  practical  problems  are  in  a  class  of  localized 
dalamination  whose  growth  behavior  is  two-dimensional  in  nature.  The  same 
energy  release  rate  concept  can  still  be  applied,  provided  a  three  dimen¬ 
sional  stress  solution  technique  is  made  available  for  the  calculation  of 
the  energy  release  rate  along  the  contoured  delamination  front.  This  topic 
is  still  new;  future  research  is  required. 
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4,  Experimental  results  from  fatigue  tests  seem  to  suggest  that  mode-I 
delamination  is  generally  static  load  related  failure,  while  mixed-mode 
delamination  is  fatigue  load  related  failure.  Correlations  between  test 
results  and  analysis  also  reveal  this  tendency.  The  fundamental  mechanisms 
of  fatigue  failure  remain  yet  to  be  investigated. 

5.  On  the  basis  of  a  limitted  fatigue  delamination  data,  a  simple  power 
law  for  the  growth  rate  was  applied  in  several  case  studies.  Though  still 
tentative  in  nature,  the  growth  law  seemed  to  show  a  considerable  generality 
for  fatigue  growth  simulation.  Again,  more  research  is  neened  here. 
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Figure  4,  Rationale  For  Predicting  the  Onset  Load  of  Delamination 
According  to  Fracture  Mechanics  Criterion. 


Figure  5.  Experimental  Edge  Delamination  Growth  Curves  Under  Uniaxial  Tension.  Ref.  £30]. 


NADC-84017-60 


Release  Rate  Coefficient  Function  C  for  Delami nation  in  C±25/90D  Laminate 


on  C_  for  Delamination  in  £±25/90 []  Laminate 


Applied  Laminate  Stress  3^  (compression) 

Figure  8.  Experimental  Edge  Delamination  Growth  CurvesUnder  Uniaxial  Compression.  Ref.  C3l* 


Figure  9,  Through-Thickness  Distribution  for  the  Interlaminar  Stress  Near  Free 


Figure  10.  Through-Thi ckness  Distribution  for  the  Interlaminar  Stress  d  Near  Free 


Through-Thickness  Distribution  for  the  Interlaminar  Stress  r  Neax  Free 


Through-Thickness  Distribution  for  the  Interlaminar  Stress  r  Near  Free 


Through-Thickness  Distribution  for  the  Interlaminar  Stress  r  Near  Free 


Release  Rate  Coefficient  Function  for  Delamination  Under  Compression 


Release  Rate  Coefficient  Function  C  for  del ami nation  Under  Compression. 


log(cy.  to  failure) 

Figure  19.  Tension  And  Torsional  Fatigue  Test  Results  of  3501-06  Epoxy  Resin.  Ref. [32] 
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Experiment  and  Predicted  Edge  Delamination  Growth  Under  Tension  Fatigue.  Ref.  C30] 
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Figure  21.  Comparison  Between  the  Experimental  and  the  Predicted  Delamination  Growth  Results 


compression  fatigue 


Figure  22.  Comparison  Between  the  Experimental  and  the  Predicted  Delamination  Growth  Results. 
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25*  Stresses( non- trivial)  in  the  45  Layer  Immediately  Above  the  Laminate  Mid-plane 


Figure  26.  (a)  Nodal  Release  Sequence  Near  Interface  Corner} 

(b)  Computed  Cfl  Coefficient  According  to  Prescribed 

Nodal  Release  Sequence. 
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